Convergent Calculation of the Asymptotic Dimension of Diffusion Limited 
Aggregates: Scaling and Renormalization of Small Clusters 
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Diffusion Limited Aggregation (DLA) is a model of fractal 
growth that had attained a paradigmatic status due to its sim- 
plicity and its underlying role for a variety of pattern forming 
processes. We present a convergent calculation of the frac- 
tal dimension D of DLA based on a renormalization scheme 
for the first Laurent coefficient of the conformal map from 
the unit circle to the expanding boundary of the fractal clus- 
ter. The theory is applicable from very small (2-3 particles) 
to asymptotically large (n — * oo) clusters. The computed 
dimension is D = 1.713 ± 0.003. 

The model of Diffusion Limited Aggregation was pre- 
sented in 1981 by Witten and Sander 0] as a computer 
algorithm. The process begins with fixing one particle at 
the center of coordinates in <i-dimcnsions, and follows 
the creation of a cluster by releasing random walkers 
from infinity, allowing them to walk around until they 
hit any particle belonging to the cluster. The growing 
cluster appears to be a random fractal. The model has 
attracted enormous attention as an example for the spon- 
taneous creation of fractal objects in nature, and also as 
a paradigm for a family of related "harmonic" problems 
that have to do with the solution of the Laplace equa- 
tion V 2 p = with appropriate boundary conditions on 
moving, ramified boundaries. Among such problems are 
dielectric breakdown flows in porous media ( with 
D'arcy's law v = — Vp and V • v = 0) ||, electrochemi- 
cal deposition Q etc. 

Numerical estimates of the fractal dimension D of DLA 
H turned out to converge extremely slowly with the num- 
ber of particles n of the cluster, leading even to specu- 
lations @ that asymptotically the clusters were plane 
filling (i.e. D = 2 in two dimensions). To date there 
are still no controlled calculations of the fractal dimen- 
sion of DLA. The aim of this Report is to close this gap 
for DLA in two dimensions. We propose a renormaliza- 
tion procedure that culminates with an integral equation 
whose solution determines the dimension D of DLA. 

To derive the wanted integral equation we use the con- 
formal theory that was developed recently [f^-jlTf for frac- 
tal growth patterns. In this theory one considers the 
conformal map from the exterior of the unit circle in the 
complex plane to the exterior of the growing fractal clus- 
ter. At the basis of this approach lies the understanding 
that once a fractal object is well developed, it is extremely 
difficult to find a conformal map from a smooth region to 
its boundary, simply because the conformal map is terri- 
bly singular on the tips of a fractal shape. The derivative 
of the inverse map is the growth probability for a random 



walker to hit the interface (known as the "harmonic mea- 
sure" ) which has been shown to be a multifractal measure 
p2| characterized by infinitely many exponents |D|jl4j . 
Accordingly, in the present approach one grows the clus- 
ter by iteratively constructing the conformal map start- 
ing from a smooth initial interface. Consider ^^(w) 
which conformally maps the exterior of the unit circle 
e %e in the mathematical iw-plane onto the complement of 
the (simply-connected) cluster of n particles in the phys- 
ical z-plane J7|-pl|. The unit circle is mapped onto the 
boundary of the cluster. The map $(™)(it;) is made from 
compositions of elementary maps (f>\ $, 



$<»>(«,) =$ (n - 13 (K,A») 



(1) 



where the elementary map 4>\. g transforms the unit circle 
to a circle with a "bump" of linear size vA around the 
point w = e l6 . An example of a good elementary map 
4>\,e was proposed in Q, endowed with a parameter a in 
the range < a < 1, determining the shape of the bump. 
We employ a — 1/2 which is consistent with semicircular 
bumps. Accordingly the map <f>( n \w) adds on a new 
bump to the image of the unit circle under $( n-1 )(u;). 
The bumps in the z-plane simulate the accreted particles 
in the physical space formulation of the growth process. 
Since we want to have fixed size bumps in the physical 
space, say of fixed area Aq, we choose in the nth step 



A n — 



A„ 



(2) 



The recursive dynamics can be represented as iterations 
of the map 0a„,0„ {w ), 



$(™)( W ) = 0Ai,0i ° <f>\ 2 ,e 2 ° • ■ ■ ° 0A„,0„M 



(3) 



The difference between various growth models will man- 
ifest itself in the different itineraries {9\ ■ ■ ■ 6 n } ||lo|j . To 
grow a DLA we have to choose random positions n . This 
way we accrete fixed size bumps in the physical plane ac- 
cording to the harmonic measure (which is transformed 
into a uniform measure by the analytic inverse of $(")). 
The DLA cluster is fully determined by the stochastic 
itinerary {&k}k=i- I n Fig. ?? we present a typical DLA 
cluster grown by this method to size n =100 000. 

This method affords us analytic power that is lacking in 
the original computer algorithm jl|; The conformal map 
<j>( n ) ( w ) i s represented in terms of its Laurent expansion: 



$ (n) H = f[ k) w - 



rliW 



jp( n ) -2 



(4) 
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The recursion equations for the Laurent coefficients of 
$(")(«;) can be obtained analytically, and in particular 
one shows that 

n 

F[ n) = J| y/l + A fc , (with the choice a = 1/2) . (5) 



fe=i 



The first Laurent coefficient has a distinguished role 
in determining the fractal dimension of the cluster, be- 
ing identical to the Laplace radius which is the radius 
of a charged disk having the same field far away as the 
charged cluster S. Moreover, defining R n as the mini- 
mal radius of all circles in z that contain the n-cluster, 
one can prove that p5[ 



Rn < 4F n 



(n) 



(6) 



Of course, for every realization {Oi}f =1 the first Laurent 
coefficient is a random number depending on and 

on Ao- ft is thus natural to consider the mean of F^ 
over all the possible realization of growth: 



(n) 



(ff>>(Ao)=/ d6 x --- d0 n F^\{0^ =1 ,X o ) (7) 



In light of Eq. (g) one expects that for sufficiently large 
clusters 



(F 1 (n) )(A )~ V^n 1 ' , 7^00 



(8) 



But this is true only for very large values of n. For arbi- 
trary values of n we offer the following proposition, which 
is central to our developments: 

proposition : For Ao of 0(1) {F[ )(Ao) is a scaling 
function of the single variable x = VAo (n + a) 1 ' ' D , where 
n takes on integer values n = 1, 2, 3 • • • and a is a weak 
function of Ao, taking on values of the order of unity. 

The origin of the parameter a is evident - it stems 
from the fact that for varying values of Ao the unit circle 
around which we grow the bumps is contributing to the 
Laplace radius, changing the effective value of n. We will 
show however that to very good approximation a can be 
taken as constant. 

We demonstrate the proposition by a direct calculation 

of (F^HAo). In Fig.2 panel a we show {F[ n) )(\ Q ) for 
ten different values of Ao as a function of n. In panel b we 
show the same data as a function of the scaling variable 
x using a = 0.2, D = 1.71. The data collapse is evident. 

We stress three points: (i) The scaling function appears 
to exist for all values of n starting from n = 1. It is 
a non-linear function of the scaling variable, with the 
attainment of the linear regime (JsJ) not in sight, (ii) a was 
taken as a Ao-independent constant; Even better data 
collapse could be obtained with a Ao-dependent a, and 
we return to this issue below, (iii) The attainment of 
data collapse requires a value for D; we have used D = 
1.71 but close-by values would have done equally well 



by changing a a bit. Thus we cannot propose the data 
collapse as an accurate method of computing the fractal 
dimension D. For this purpose we derive now an integral 
equation from which both a and D can be computed in 
a controlled fashion. 

Having a scaling function in mind we think about a 
real-space renormalization group procedure, in which we 
change the number of particles in the cluster and their 
size such as to keep the cluster invariant. For n very large, 
when the radius of gyration is linear in -\/A~o n 1 / , we can 
change the size of the particles VAo by a factor of 2, and 
their number n by a factor of 2 D . In the asymptotic 
regime such a renormalization will leave R n invariant. 
For small values of n we make use of the discovery that 
(F\ )(Ao) is a scaling function of the single variable x = 
y/\o(n + a) 1 / . We demand that upon renormalization 
the average Laplace radius remains invariant. In other 
words, the fixed point condition can be written as the 
following integral equation: 



2-k r2-K 

d9x--- de n F[ n \\ Q ) 



n Jo 

2k r 2ir 



Mx-- - I d9 n Fl n) (X ), 
o Jo 



(9) 



where for any n > n the equation is satisfied by the 
unique value of Xq that solves the equation 



Ao — V A. 



n + a 
n + a 



l/D 



(10) 



The way to compute D is then obvious - one computes 
the integral on the LHS of Eq.{j9j) for some value of n, 
and then finds the unique value of A for which the RHS 
with n > n equals the LHS. Then 



D = 2 



ln(n + a) — ln(n + a) 
In Aq — In Ao 



(11) 



One should stress that Eq.(g) is all explicit - there is 
no simulation or randomness left at this point - the inte- 
grand on both sides is an explicit function of {6*i, • • • , 8 n } 
or {6»i, • • • , 8n) through Eqs.(|) and (§). We note that 
this method of calculation is fundamentally different from 
the standard method of log-log plots of the radius of the 
cluster vs. n [^|j6). These rely on the proportionality of 
R n and n x l D . Here we may use small values of n as we 
depend neither on the asymptotic linearity of (|8|) , nor on 
self-averaging. 

Clearly the integral equation will be useful for an ac- 
tual determination of D only if it converges to the fixed 
point quickly upon increasing n. Otherwise the calcula- 
tion of the multidimensional integral will become very 
cumbersome, maybe even more time consuming than 
standard numerical simulations. In the rest of this Re- 
port we demonstrate that the calculation actually con- 
verges very quickly and present the determination of D. 
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Firstly we determine the value of a. We use Ao = 1, 
and solve Eq.(||) for the simplest case n = 1, n = 2, 
demanding 

(f[ 2) )(\ = i) = V2(VT+T 2 ) = f\ (1) (a ) = v/TTa^ (: 

{2 + a) 1 / = (l + a) 1 / ^ (: 

Computing explicitly (taking D = 1.71) we found 

+ A 2 ) w 1.241 from which followed a ps 0.146. To 
bracket the calculations we have employed this value to- 
gether with a = and a = 1. 

In solving Eq.(|J) we chose invariably n = n + 1, and 
Ao = 1. The results are shown in Fig. 3, in the form of the 
value of D vs. (n+a)^ 1 . Fitting the best nonlinear curve 
to the data with a — 0.146 we find that it extrapolates 
for (n + a)^ 1 — ► to D = 1.7150.... Obviously the calcu- 
lation with a — and a — 1 bracket this from above and 
from below. Nevertheless the best nonlinear fits to these 
data extrapolate to very close values for (n + a)^ 1 — > 
(see Fig. 3). Taking all the data together we can present 
a final value of D = 1.713 ± 0.003. We cannot overstress 
the fact that these results were obtained from solving 
Eq.(||) with n values ranging from n = 1 to n = 20. We 
believe that this represents a major advance compared 
to traditional estimates of the dimension of DLA. 

To improve the results even further, and to remove the 
curvature in the line of D vs. (n + a)^ 1 we can endow a 
with a weak Ao dependence. Using the function 

a(A ) = 2.04811 - 4.08071 VA^ + 2.23446A , (14) 

and using the appropriate value of a(Ao) in Eq.(|i"T|) we 
obtain solutions with essentialy constant D for all values 
of n > 2. The values of D computed are shown in Fig. 4 
as a function of n. With this data we can state without 
going to the limit n ->• oo that D = 1.713 ± 0.005. This 
result is invariant with respect to changing n in Eq.(^|). 
It can be computed equally well from n = 3 and h — 4 
as from and n = 19 and n = 20. We should of course 
stress at this point that the function a(Ao) was not com- 
puted from first principles; we hope that further theoret- 
ical progress will shed light on how to achieve a-priori 
determination of the functional form a(Ao). 

The most important questions that we need to address 
now are: (i) why the classical numerical estimates |^|,^) of 
the fractal dimension of DLA converge so slowly, whereas 
here we can get an excellent estimate of D even with n 
of the order of unity, (ii) Is the collapse of the scal- 
ing data for small n exact or an excellent approxima- 
tion. The answer to (i) is that in standard numerical 
experiments the radius of gyration of the grown cluster 
was plotted in log-log coordinates against the number of 
particles, with D estimated from the slope. Examining 
our scaling function (Fi)(x) (see Fig. 2) we note the 
slow crossover to linear behavior, which may not be fully 
achieved even for extremely high values of n. In this re- 
spect we understand that a reliable estimate of D from 
radius of gyration calculation requires inhuman effort, as 



was indeed experienced by workers in the field [pi . In the 
present formulation the appearance of the asymptotic D 
as a renormalization exponent already at early stages of 
the growth allows a convergent calculation. The answer 
to (ii) has to await a first principle theory for a(Ao). 
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FIG. 1. Panel a: the mean (F^) as a function of n for 
different values of A . The data shown pertain to the 10 
values of Ao from 0.6 (below) to 1.5 (above). Panel b: the 
same objects plotted as a function of the scaling variable x 
with a = 0.2 and D = 1.71. 
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FIG. 2. The estimated value of the fractal dimension D 
computed from the solution of the integral equation (9) . The 
three sets of points pertain to a — (upper), a = 0.146 
(middle) and a = 1 (lower). The lines are the best quadratic 
fits to the data. The extrapolated values at l/(n + a) — 
are 1.7158.., 1.7150.. and 1.7111.. respectively, leading to our 
final estimate D = 1.713 ± 0.003 
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FIG. 3. The estimated value of the fractal dimension D 
computed from the solution of the integral equation (9), using 
the function a(A ) of Eq.(14). Without going to the limit 
n — > co we can state D = 1.713 ± 0.005 
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This figure "deterDLA100000.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0008053vl 



